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A  SIMPLE  METHOD  OF  DETERMINING  STRESS  IN 
CURVED  FLEXURAL  MEMBERS* 


I.  Introduction 

1.  Calculation  of  Stresses  in  Curved  Beams. — Curved  members  sub- 
jected to  bending,  such  as  crane  hooks,  frames  for  punches  and  drill 
presses,  members  in  railway  trucks,  chain  links,  etc.,  are  important  ele- 
ments in  many  engineering  machines  and  structures.  It  is  generally 
understood  that  the  ordinary  flexure  formula,  which  gives  reliable  values- 
for  the  primary  stresses  in  straight  beams  or  flexural  members,  gives 
values  for  stress  that  may  be  much  too  small  when  applied  to  curved 
flexural  members.  Nevertheless,  this  straight-beam  formula  is  still  widely 
used  in  the  design  of  curved  flexural  members. 

The  main  reasons  for  applying  the  straight-beam  formula  to  curved 
beams  are,  probably,  that  the  curved-beam  formula  is  much  more  com- 
plicated and  difficult  to  apply  than  is  the  straight-beam  formula  and 
that,  in  the  case  of  cast  frames,  the  deviation  from  the  assumed  con- 
dition of  homogeneity  of  material,  as  well  as  the  presence  of  initial 
stresses,  causes  uncertainty  as  to  the  significance  of  the  calculated 
stresses.     These  conditions  will  be  discussed  briefly  later  in  this  circular. 

Due  to  the  shorter  length  of  the  fibers  on  the  concave  (or  inner) 
side  of  a  curved  beam  the  unit  strain  (and  hence  also  the  unit  stress 
within  the  proportional  limit)  is  greater  than  that  given  by  the  straight- 
beam  formula ;  and  the  stress  on  the  convex  (or  outer)  side  is  less.  Like- 
wise the  neutral  axis  does  not  pass  through  the  centroid  of  the  cross- 
section  of  the  beam  but  lies  nearer  to  the  concave  or  inner  side. 

This  increase  in  stress  due  to  the  curvature  of  the  beam  may  be  of 
great  importance,  especially  in  members  in  which  the  stress  is  repeated 
many  times,  as  is  frequently  the  case. 

Winkler f  was  perhaps  the  first  engineer  to  point  out  clearly  the 
difference  between  the  stress  distribution  in  a  straight  beam  and  in  a 
curved  beam,  and  to  emphasize  the  engineering  importance  of  the  prob- 
lem. Bach|  extended  the  analysis  and  also  obtained  valuable  experi- 
mental results.  The  formula  resulting  from  this  analysis,  which  will  be 
referred  to  as  the  Winkler-Bach  formula,  §  has  gained  rather  wide  accep- 
tance; it  may  be  written 

♦The  discussion  herein  presented  applies  only  to  members  subjected  to  symmetrical  loading,  that  is, 
where  the  plane  of  the  loads  contains  a  principal  axis  of  inertia  of  each  cross-section,  and  curved  only  in 
the  plane  of  loading. 

tFor  Winkler's  memoir  see,  "Formanderung  und  Festigkeit  Gekrummter  Korper  ins  besondere  der 
Ringe."  Der  Civilingenieur  Bd.  IV  S.  232-246. 

tSee  C.  Bach,  "Elastizitiit  und  Festigkeit."  See  also  "The  Strength  of  Chain  Links,"  Univ.  of 
111.  Eng.  Exp.  Sta.  Bui.  18. 

§This  formula  is  derived  in  Appendix  C. 
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for  combined  direct  stress  and  bending. 
S    =  unit  stress  at  fiber  in  question 

M  —  bending  moment  (positive  when  it  increases  curvature,  neg- 
ative when  it  decreases  curvature) 
=  area  of  cross-section 

=  distance  from  center  of  curvature  to  centroidal  axis  of  sec- 
tion 
=  distance  from  centroidal  axis  to  fiber  in  question  (negative 
when  measured  toward,  and  positive  when  measured 
away  from  center  of  curvature) ;  for  the  extreme  fiber  the 
value  of  y  is  denoted  by  c  as  in  Fig.  1. 
=  a  property  of  the  cross-section  defined  by  the  expression, 
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Si  =  axial  or  direct  unit  stress 

This  formula  has  been  shown  to  give  results  that  agree  well  with  the 
stresses  found  from  experimental  strain  measurements  of  curved  beams 
made  of  steel.*  The  equation  a'so  satisfies  the  limiting  condition  of 
reducing  to  the  ordinary  flexure  formula  when  the  radius  of  curvature 
becomes  very  large;  that  is,  when  the  beam  becomes  practically  straight. 

The  formula,  however,  although  satisfactory  theoretically,  has  a 
serious  disadvantage  because  of  its  difficulty  of  application,  the  difficulty 
being  in  evaluating  the  property  Z  of  the  cross-section.  The  evaluation 
of  this  property  of  the  area  by  an  analytical  method  is  complicated  for 


Ax/'s  of  C£sry&/-£s/-<£' 

Any         trji* 
F/frer-    r 


/\x/s 


Cen/ro/a'a'/ 
Ax/s 


Secf/or?  A-B 

Fig.  1.   A  Typical  Curved  Flexukal  Member 


*"Tests  and  Theory  of  Curved  Beams,"  Window  and  Edmonds,  A.  S.  M.  E.  paper,  1926. 
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all  sections  except  the  circle  and  the  rectangle;  and  even  for  the 
simpler  ordinary  composite  sections,  as  for  instance  the  I  or  T  sec- 
tions, the  expression  becomes  too  complicated  for  quick  or  accurate 
solution.* 

A  possible  solution  of  this  difficulty  of  evaluating  Z  is  to  use  graph- 
ical methods.  One  graphical  method  is  the  "link  rigidity"  method,  f  and 
another,  which  the  writers  have  found  to  be  more  accurate,  is  the  "line 
polygon"  method. %  Either  method  is  long  and  susceptible  to  those 
errors  inherent  in  graphical  work,  which  in  this  case  are  serious,  since 
Z  is  a  small  quantity  for  any  section,  rarely  exceeding  0.8  for  a  very 
sharp  curvature,  and  becoming  less  than  0.005  for  a  relatively  large  curva- 
ture. Experience  with  graphical  methods  of  determining  Z  has  convinced 
the  writers  that  the  error  that  may  be  expected  is  at  least  ten  per  cent, 
and  may  be  considerably  greater.  §  Graphical  methods  of  determining 
Z  may  be  considered  as  a  last  resort  when  analytical  expressions  cannot 
be  obtained,  or  are  very  complicated. 

A  way  of  evading  these  difficulties  has  frequently  been  suggested; 
namely,  to  use  the  straight-beam  flexure  formula  with  a  correction 
factor.    The  expression  for  pure  bending  in  a  curved  beam  would  then  be 

S  =  K  (Mc/I)  (3) 

and  for  combined  direct  stress  and  bending 

3  =  S,  +  K  (Mc/I)  (4) 


where  K  = 


M 

"i  1 1(  c  VI 

aR 

L     '  z\R  +  c)_ 

Mc 
I 

(5) 


No  values  of  K  seem  to  be  available  except  for  circular  and  rectan- 
gular sections,  and  although  it  has  generally  been  assumed  that  K  is 
approximately  the  same  for  sections  of  all  shapes,  no  careful  study  of 
the  value  of  K  seems  to  have  been  made. 

The  object  of  this  investigation,  therefore,  was  to  determine  accur- 
ately values  of  the  correction  factor  K  for  variously  shaped  sections  and 
for  degrees  of  curvature  commonly  found  in  curved  flexural  members; 
and  thus  to  make  available  the  results  of  the  curved-beam  theory  without 
involving  the  usual  difficulties  and  errors  encountered  in  determining 

*See  Appendix  A  for  analytical  expressions  for  Z. 
tFor  description  see  "Strength  of  Materials,"  A.  Morley,  -p.  383. 
JSee  Appendix  B. 

§For  a  discussion  of  the  errors  likely  to  occur  in  the  value  of  Z  see  "Bending  Stresses  in  Hooks  and 
Other  Curved  Beams,"  by  A.  Morley,  Engineering  (London)  Vol.  98,  Sept.  11  and  25.  1914. 
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the  value  of  Z  and  in  using  the  more  involved  curved-beam  formula. 
The  location  of  the  neutral  axis  in  a  curved  beam,  a  knowledge  of  which 
is  sometimes  of  real  value,  is  also  given  for  various  cross-sections  and 
degrees  of  curvature. 

2.  Acknowledgments. — The  work  herein  presented  was  originally 
undertaken  as  a  part  of  a  graduate  course  at  the  University  of  Illinois 
under  Fred  B.  Seely,  Professor  of  Theoretical  and  Applied  Mechanics, 
and  was  completed  under  his  direction.  The  writers  gratefully  acknowl- 
edge his  many  helpful  suggestions  and  generous  assistance. 

Credit  is  due  Mr.  Albert  E.  Hershey,  Research  Assistant  in  Mechan- 
ical Engineering,  for  the  graphical  solution  of  Z  and  the  analytical 
expressions  for  Z  given  in  the  appendices.  This  work  was  done  previously 
as  part  of  a  graduate  course  under  Professor  Seely 

II.  Correction  Factor  for  Use  With  Ordinary 
Flexure  Formula 

3.  Values  of  Correction  Factor. — The  value  of  the  correction 
factor  K  has  been  computed  for  the  inside  fiber  A,  and  the  outside  fiber 
B  (Fig.  1)  for  12  of  the  more  common  sections,  and  for  10  different 
degrees  of  curvature  in  each  case.  Every  effort  has  been  made  to  obtain 
accurate  results.  The  values  of  Z  were  obtained  from  the  analytical 
expression  by  use  of  a  calculating  machine,  and  in  all  cases  were  checked 
by  means  of  graphical  methods.     The  resulting  values  of  K  are  given 

R 

in  Table  1.    \  allies  of  K  for  intermediate  values  of   -  may  be  found 

without  appreciable  error,  by  ordinary  linear  interpolation. 

The  range  of  the  degree  of  curvature  has  been  so  chosen  as  to 
include  the  two  extreme  conditions.  In  the  one  case  the  curvature  is  so 
sharp  as  to  approach  the  condition  of  localized  stress,  under  which  the 
mass  action  of  the  material  and  the  distribution  of  stress  assumed  n 
the  ordinary  formulas  of  mechanics  no  longer  hold  ;*  and  in  the  other  case 
the  curvature  is  so  slight  that  the  uncorrected  flexure  formula  may  be 
applied  without  serious  error. 

It  is  evident  that  there  is  considerable  variation  in  the  value  of  K 

for  different  sections  for  values  of  —  of  2.0  or  less;  that  is,  for  beams 

c 

in  which  the  curvature  is  of  real  importance. 

Table  1  gives  values  of  K  for  sections  that  may  be  considered  as 

typical.    Fortunately  there  are  sections  for  which  K  is  independent  of 

the  breadth  of  the  section,  as  noted  in  Table  1.    For  other  sections  K 

*See  F.  B.  Seely,  "Resistance  of  Materials,"  Appendix  III. 
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Table  1 
Values  of  K  for  Different  Sections  and  Different  Radii  of  Curvature 
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Table  1  (Continued) 
Values  of  K  for  Different  Sections  and  Different  Radii  of  Curvature 
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l/„  as  o''sfa'/?<?e   fro/77    cer>f/-0/a'cr/  orgy's  So  /7ec/fr&/  o'^-z's. 
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Table  1  (Concluded) 
Values  of  K  for  Different  Sections  and  Different  Radii  of  Curvature 


Secf/or? 

/? 

Fcrcfor  /V 

* 

C" 

/r)s/a/e  F//?e/~ 

Ou/s/ae  F/ber 

/.2 
/.4 

3.Z8 
2.3/ 

0.58 
0.64 

0.266/? 
0.182/? 

I 

i 

/.6 

1.89 

0.68 

0.134  F 

/.8 

1.70 

0.7/ 

0./04/? 

W^&-jOyy. 

1 

2.0 

f.5? 

0.73 

0.083/? 

| 

3.0 

1.3/ 

0.8/ 

0.038/? 

4.0 

1.2/ 

0.86 

0.020/? 

I 

6.0 

U3 

0.90 

0.0087/? 

U      /? 

_J 

8.0 

/JO 

0.92 

0.0049/? 

// 

/O.O 

/.07 

0.93 

0.003//? 

r — -*4t — i 

1 

/.z 

/.4 

/.6 

2.63 
/.97 
/.66 

0.68 
0.73 
0.76 

0.399/? 
O.280/? 
0.205/? 

///y// 

1 

/.8 

/.5/ 

0.78 

0/S9/? 

>  f+ 

— <?/•— 

— f    ! 

Z.O 

1.43 

0.80 

0./Z7/? 

,        i 

3.0 
4.0 
6.0 

1.23 
US 
f.09 

0.86 
0.89 
0.92 

0.058/? 
0.03//? 
0.0/4/? 

1 

~c-A 

See  a/so  I  secf/os?  Ta£/e  2 

8.0 

i.0? 

0.94 

O.0O76/? 

/2 

/O.O 

/.OS 

0.95 

0.0048/? 

*y„  is  o//s/a'/?ce  fro/7?  ce/?fs-o/a/a'/  avr/'s  fo  /7£>£/f/-a'/  <?/c/s. 


varies  appreciably  for  varying  section  proportions  as  indicated  by  the 
values  for  the  trapezoidal,  T,  and  I  sections;  but  it  is  felt  that  the  value 
of  the  correction  factor  for  any  section  likely  to  be  used  in  a  curved 
flexural  member  may  be  estimated  with  sufficient  accuracy  directly  from 
Table  1. 

4.  Empirical  Formula  for  K. — Although  Table  1  may  be  cons  dered 
as  sufficiently  indicating  the  approximate  value  of  a  correction  factor  for 
any  section  used  in  curved  members,  nevertheless  it  seemed  desirable 
to  obtain  an  empirical  formula  for  K  on  the  basis  of  the  values  in  Table 
1,  and  to  check  the  values  obtained  by  the  use  of  this  formula  when 
applied  to  sections  different  from  those  in  the  table  against  values 
calculated  in  the  same  way  as  those  in  the  table.  The  formula  developed 
is  as  follows: 


K  =  1.00  +  0.5 


be2  \-R 


1 


c  +  R 


(6) 


(use  1.05  instead  of  0.5  for  circular  and  elliptical  sections) 
/    =  moment  of  inertia  as  used  in  straight-beam  formula 
b     —  maximum  breadth  of  section 

c     =  distance  from  centroidal  axis  to  inside  fiber,  i.e.,  to  the  ex- 
treme fiber  nearest  the  center  of  curvature 
R    =  radius  of  curvature  of  centroidal  axis  of  beam 
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It  will  be  observed  that  the  value  of  K  in  this  formula  depends  on  the 
distribution  of  the  area  as  involved  in  the  expression  f  —  ),  and  also  on 

r  i    ,n 

the  ratio     R  —  c       R     wnicn  takes  account  of  the  curvature  of  the 

section. 

In  spite  of  its  simplicity  this  formula  for  K  gives  very  satisfactory 
results  for  all  sections  (except  the  triangle)  in  Table  1  throughout  the 
entire  range  of  curvatures.  Values  of  K  obtained  by  use  of  this  formula 
will  agree  with  the  computed  values  for  any  section  in  Table  1,  with  the 
exception  of  the  triangle,  with  a  maximum  error  of  5  per  cent,  the 
average  error  being  considerably  less  than  this.  In  the  case  of  the  tri- 
angle the  error  is  as  great  as  10  per  cent  in  one  case,  but  it  is  a  section 
rarely  used,  and  was  included  in  the  table  only  as  an  extreme  case  of  the 
trapezoid.  Hence  no  particular  effort  was  made  to  fit  its  K  curve  with 
the  formula. 

This  empirical  formula  has  been  tested  for  numerous  section  pro- 
portions not  included  in  Table  1,  and  a  maximum  error  of  5  per  cent 
seems  to  cover  all  cases.  An  error  of  5  or  even  10  per  cent,  however,  is 
distinctly  less  than  that  with  which  the  load  on  a  member  can  be  pre- 
dicted, or  the  assumption  of  homogeneity  of  material  assured.  Further, 
the  error  in  calculating  Z  by  the  regular  method  can  easily  exceed  10 
per  cent  unless  five-place  accuracy  is  used  in  its  determination,  and  even 
then  the  lengthy  computation  involved  makes  mistakes  readily  possible. 

The  formula  as  given  applies  to  the  inside  fiber  only,  i.e.,  to  the 
extreme  fiber  on  the  concave  side.  Since  the  stress  in  the  extreme  fiber 
on  the  convex  side  is  always  less  than  that  given  by  the  straight-beam 
flexure  formula,  and  since  the  significant  stress  is  nearly  always  that 
at  the  inside  fiber,  no  attempt  has  been  made  to  derive  an  empirical 
formula  from  which  the  K  for  the  outside  fiber  can  be  obtained. 

It  will  be  observed  that  the  empirical  formula  satisfies  the  condition 
that  K  approaches  unity  as  the  radius  of  curvature  becomes  very 
large,  and  the  member  approaches  a  straight  beam.  Likewise,  as  the 
radius  of  curvature  becomes  small  (approaches  zero)  the  value  of  K  be- 
comes indefinitely  large  and  hence  the  calculated  stress  at  the  inside 
fiber  becomes  infinitely  large,  which  is  consistent  with  the  Winkler-Bach 
formula. 

The  values  of  K  for  some  of  the  sections  listed  in  Table  1  are  shown 
graphically  in  Fig.  2.  The  variation  in  the  value  of  K  corresponding 
to  the  sharper  curvatures  is  evident. 

5.  Limitations  and  Assumptions. — It  should  be  noted  that  the 
values  of  correction  factors  in  Table  1  should  be  used  in  calculating  the 
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Fig.  2.   Variation  of  K  with  —  for  Three  Sections 

c 


stress  only  in  the  extreme  fibers  of  the  curved  beam.  However,  this  is 
where  the  maximum  stress  occurs.  In  order  to  show  the  distribution  of 
stress  over  the  cross-section,  and  to  illustrate  the  fact  that  the  ordinary 
flexure  formula  with  these  correction  factors  applies  only  to  the  extreme 
fiber,  the  sketch  in  Fig.  3  is  shown. 

One  of  the  fundamental  assumptions  on  which  all  mathematical 
determination  of  stresses  in  members  is  based  is  that  the  material  is 
homogeneous  and  isotropic,  and  hence  that  there  are  no  discontinuities 
in  the  stress  distribution  on  the  cross-section.  But,  although  no  material 
satisfies  this  assumption,  the  values  given  by  the  usual  design  formulas 
are  in  general  satisfactory  at  least  for  ductile  material  subjected  to 
static  loads.  In  other  words,  even  though  the  calculated  stress  is  only 
an  average  or  statistical  value  of  the  stress  in  the  ne'ghborhood  of  the 
point  at  which  a  knowledge  of  the  stress  is  desired,  it  is  nevertheless  the 
significant  stress  for  most  materials  subjected  to  static  loads.  Thus  the 
Winkler-Bach  curved-beam  formula  gives  the  significant  stress  for  most 
materials.  The  question  is  sometimes  raised,  however,  as  to  whether  the 
value  calculated  from  the  curved-beam  formula  is  more  significant  in  the 
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P  P7c 


5= 


P   P7<r 


Fig.  3.    Stress  Distribution  in  a  Curved  Flexural  Member 
as  Given  by  Three  Formulas 


case  of  cast  material,  particularly  cast  iron,  than  is  the  value  given  by 
the  straight-beam  formula.  In  the  calculation  of  the  stress  in  curved 
beams  by  means  of  the  uncorrected  straight-beam  formula,  the  so-called 
factor  of  safety  must  cover  both  the  uncertainty  of  the  material  and  the 
error  arising  due  to  the  use  of  an  incorrect  ormula.  (The  stress  in  the 
curved  beam  may  be  more  than  3  times  as  great  as  the  value  given  by  the 
straight-beam  formula,  see  Table  1).  This  method  of  design  was  perhaps 
justified  because  of  the  difficulties  involved  in  the  use  of  the  Winkler- 
Bach  formula.  The  use  of  the  straight-beam  formula  with  a  correction 
factor  which  is  easily  obtained  would  seem  to  make  available  for  practi- 
cal use  the  advantages  of  the  curved-beam  theory,  and  hence  to  give  a 
more  rational  basis  for  design.  Further,  this  procedure  should  help  to 
reduce  the  number  of  failures  of  curved  flexural  members,  particularly 
when  such  members  are  subjected  to  impact  and  to  repeated  stress. 

The  question  is  also  raised  as  to  whether  the  curved-beam  theory 
should  be  applied  to  a  curved  member  at  a  section  where  the  inner  and 
outer  surfaces  are  not  concentric  curves.  The  curved-beam  formula  is 
based  on  the  same  assumption  in  this  respect  as  is  the  straight-beam  for- 
mula; namely,  that  the  beam  has  a  constant  cross-section;  or,  at  least, 
that  the  change  in  cross-section  along  the  beam  is  very  gradual.  At  the 
corner  sections  of  some  frames  where  nearly  a  right  angle  turn  is  made, 
both  the  straight-beam  and  the  curved-beam  formulas  would  be,  per- 


DETERMINING  STRESS  IN  CURVED  FLEXURAL  MEMBERS  25 

haps,  considerably  in  error,  but  in  most  cases  the  value  of  the  stress  given 
by  the  curved-beam  formula  would  be  the  better  indication  of  the 
significant  stress. 

III.  Conclusions 

6.  Summary. — The  following  brief  summary  is  offered: 

(1)  The  straight-beam  flexure  formu' a  is  radically  in  error  when 
appl'ed  to  curved  flexural  members  of  relatively  sharp  curvature,  and 
this  error  is  probably  responsible  for  many  failures  of  curved  flexural 
members  particularly  when  the  members  are  subjected  to  impact  and  to 
repeated  stress. 

(2)  The  Winkler-Bach  formula  for  curved  beams,  although  re- 
liable, is  difficult  of  application. 

(3)  The  stress  in  curved  beams  can  be  satisfactorily  calculated 
by  the  ordinary  straight-beam  flexure  formula  used  with  a  correction 
factor.  Values  of  the  correction  factor  for  any  sections  likely  to  be  used 
in  structural  and  machine  members  may  easily  be  obtained  by  use  of 
Table  1,  or  equation  (6). 

(4)  The  stress  computed  by  means  of  this  corrected  straight- 
beam  formula  applies  only  to  the  extreme  fibers  of  the  member.  But, 
since  the  significant  stress  usually  occurs  at  an  extreme  fiber,  the  method 
offers  a  simple,  satisfactory,  and  practical  means  of  designing  curved 
flexural  members. 

APPENDIX  A 
Expressions  for  Section  Property  Z 

As  pointed  out  in  the  d'scussion  of  the  Winkler-Bach  formula,  the 

i  C  y 

evaluation  of  the  integral  Z  =  —  -    I  ttt  —  da  is  the  chief  difficulty  en- 

aj  R  +  y 

countered  in  the  use  of  the  curved-beam  formula.  Performing  this  in- 
tegration results  in  complicated  expressions  for  all  but  the  circular  and 
rectangular  sections. 

V  .  . 

By  expanding    .,\_       iloa  series  we  get  the  general  expression 


a 


a + gr  -  iff + (f   n 


Aa 


For  symmetrical  sections,  the  summation  of  the  terms  with  odd  powers 
is  zero  and  the  series  converges  rapidly.  The  first  three  terms  are  usually 
sufficient  in  that  case,  but  where  the  series  expression  is  satisfactory,  the 
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An 


Table  2 
alytical  Expressions  for  Z 


*-  J>f£)\  i/c}4+  JiUf+ 

■-3[ffJ   +J(/?J    +7[/?/    + 


z=-/+ 


'fa.m>j 


^ 


z 
z=-, 


-ft H 


Z=-/+2 


&&»*)**(j&}^J 


-->+*(§) -*H(§r-> 


z=  -/*  ^rfV^F-  Vx^fJ 


z=-/+ 


jt>C,-3,i 


DETERMINING  STRESS  IN  CURVED  FLEXURAL  MEMBERS 


17 


Table  2  (Concluded) 
Analytical  Expressions  for  Z 


R 


Z=  -/+  ■£  [b,  /<%  (ft+c,)  +(t-b,)/ogre  f/?+oJ 
+  (b-f)/oge  (R-c3)-b  /o#e  f/?-cz)] 


J 


fjspz 


1 

1P¥m\ 


Tbe   var/ue   of  Z  for  e&cf?  of  f/?ese 
fbree   secf/or?s  /??&{/  be  foe//?*/  fro/77 
fhe  express /or?  a-boi/e  bt/  /7?a-/r/>?g' :• 
b,  =  b ,    cz  =  o,,  o/7o/  c3=c^ 

Arecr=  cr  =  2f(t-b)c,  +bcz] 


W  :\  I 


mm 


U-  c,-\-  c,~ 


/r?  fbe  express/or  for  fr?e  c/re^ua'/  I  gr/i/- 
er?  o-boi/e  mcrA-e   o4  =  o,    a/?o/  b,=t,  fber> 

-b/o#e  (ff-ca)] 
Area-=  a>  =  fc,-  (b-fjc3  +  boz 


exact  evaluation  of  the  integral  is  also  comparatively  simple,  and  is,  of 
course,  more  accurate.  For  the  unsymmetrical  sections  the  series  con- 
verges too  slowly  to  be  of  use. 

In  order  that  the  expressions  for  Z  used  in  calculating  the  values  in 
Table  1  may  be  made  available  for  those  wishing  to  use  the  curved-beam 
formula,  Table  2  is  presented.    The  evaluated  integral  is  given  for  all 
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sections  commonly  found  in  curved  flexural  members.  The  expressions 
were  obtained  from  Bach's,  "Elastizitat  und  Festigkeit,"  and  from  an 
unpublished  paper  by  A.  E.  Hershey,  and  have  been  revised  to  fit  the 
notation  of  the  figures. 


APPENDIX  B 

Graphical  Determination  of  Z  by 
Line  Polygon  Method 

Reference  has  been  made  to  the  "line  polygon"  method  of  evalua- 
ting Z  graphically.  This  method,  devised  by  A.  E.  Hershey,  is  slightly 
more  accurate  than  the  "link  rigidity"  method,  and  does  not  necessitate 
the  use  of  a  planimeter  or  other  area  measuring  device.  It  is  particular- 
ly convenient  when  the  string-polygon  method  is  used  in  locating  the 
centroidal  axis  of  the  section,  and  a  brief  description  of  it  will  be  given 
here. 

In  Fig.  4,  the  section  whose  Z  is  to  be  determined  is  enclosed  by  the 
heavy-line  curve  EABF,  etc.  Divide  the  area  into  several  strips  of  equal 
width  (about  ten  strips  should  be  used  for  accurate  results).  Draw  the 
mean  ordinates  AB.  CD,  EF,  etc.,  of  the  strips.  Then  lay  off  the  mean 
ordinates  on  a  line  A'B'l  7/.  etc.  The  mean  ordinate  of  any  strip  being 
assumed  proportional  to  the  area  of  the  strip,  the  total  length  of  this 
line  then  represents  the  area  of  the  section 

If  the  centroidal  axis  is  to  be  found  graphically,  select  a  pole  J  and 
draw  the  rays  JA',  JB',  etc.  Draw  aP  parallel  to  J  A',  and  starting  at  a, 
draw  the  strings  axi,  X\Xz,  etc.,  corresponding  to  the  rays  JB',  JC,  etc., 
to  locate  x.  Then  draw  xP  parallel  to  JK',  thus  locating  P.  The  line 
MM'  through  P  perpendicular  to  the  center  line  of  the  section  is  the 
centroidal  axis,  and  M  is  the  centroid  of  the  area. 

Now  to  obtain  a  value  of  Z,  draw  a  line  RS  parallel  to  OM  and  at 
a  distance  of  unity  from  OM,  using  the  same  scale  as  used  in  drawing 
the  section.  From  the  intersections  Vi,  F2,  etc.,  of  RS  with  the  median 
lines,  draw  lines  to  the  center  of  curvature  0.  These  lines  locate  points 
Q\,  Q-2,  etc.,  where  they  intersect  the  centroidal  axis.  Now,  at  the  ends 
of  the  vectors  A'B',  B'C ',  etc.,  in  the  vector  polygon,  erect  perpendiculars 
B'b,  C'c,  etc.,  the  lengths  of  which  are  unity  to  the  scale  used  in  drawing 
the  vector  polygon.  From  the  ends  of  these  lines  lay  off  lengths  bb' ',  cc' ', 
etc.,  equal  to  MQi,  MQ2,  etc.,  respectively.  Through  each  of  the  points 
b',  c',  d',  etc.,  thus  obtained,  draw  lines  qi,  q<i,  q%,  etc.,  parallel  to  the  vector 
polygon  A'B'C'D'E'F'G'H'K'.    The  value  of  Z  is  obtained  by  dividing 
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k'  A" 

Fig.  4.    Graphical  Determination  of  Z  by  the  Line  Polygon  Method 
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the  sum  of  the  lengths  gi,  qi,  qz,  etc.,  by  the  area  of  the  section  (represented 

by  the  sum  of  the  lengths  qi,  g2',  qz ', ),  and  subtracting  unity 

from  the  quotient.    That  is 

Z-     Sg        1 

Proof:     By  definition 

Z=  -  —  /  — L  rfa 

This  may  be  written 


z.-lfL.     «   )da_  -Lfda  +  1  f 

a  J    \  y  -\-  K/  a  J  a   J 


y  +  R 


da 


=  -  1  +  -   /"- 


Now  each  strip  is  an  elementary  area  da  and  is  at  a  distance  y  from  the 
centroidal  axis.  That  is,  MV\  =  y;  also  OM  =  /?.  Then  in  triangles 
OVSVi  and  OMQi  we  have 

In  the  triangle  A'5'6,  -^r  =  4>£-  but  &'6  =  MQ*  =    ~nZ 

bb  Bb  y+R 

R 

B'b  =  1,  and  A'B'  =  da.       Therefore  bxT  =  ql  =        ,    _  da 


-/ 


K 


da 


Hence  qx  +  52  +  93  +     -    -    -    -     =  2g  -     , 
Also,  area  a  =  g/  +  g2'  +  qz   +      -  -  -  -     =  2g' 

Therefore  Z  =  -§V  -  1 

2,g 


APPENDIX  C 

Derivation  of  the  Winkler-Bach  Formula* 

Let  the  curved  beam  DOE  (Fig.  5a)  be  subjected  to  the  loads  and 
reactions  Q.  &i  **' 

Let  AB  and  A1B1  be  two  normal  sections  of  the  beam  indefinitely 
close  together  before  the  loads  are  applied.  After  the  beam  is  loaded  the 

•The  theory  here  given  is  substantially  that  laid  down  by  Bach.  "Elastizitat  und  Festigkeit," 
Sec.  54. 
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6/     4 
(if) 

Fig.  5.  Strains  in  a  Curved  Beam 


change  in  length  of  any  fiber  PPi  is  assumed  to  be  proportional  to  its 
distance  from  the  neutral  axis  and  is  represented  by  P\Pi';  the  neutral 
surface  is  represented  by  NNi.  It  is  required  to  find  the  unit  stress  on 
any  point  P  in  a  cross-section  of  a  curved  beam  in  terms  of  the  bending 
moment  and  the  dimensions  of  the  section. 

Let  Fig.  5b  represent  the  part  ABBiAi  of  Fig.  5a  enlarged.     The 
unit  stress  S  on  a  fiber  PPi  (Fig.  5b)  at  any  distance  y  from  the  centroidal 


axis  is 


S   =    Ee 


(1) 


in  which  E  is  the  modulus  of  elasticity  of  the  material  and  e  is  the  unit 
strain  of  the  fiber 

The  bending  moment  M  causes  the  plane  AiBi  to  rotate  through 
an  angle  Add  thereby  changing  the  angle  this  plane  makes  with  the  plane 
BAC  from  dd  to  {dd  +  Add) ;  the  center  of  curvature  is  changed  from 
C  to  C,  and  the  distance  of  the  centroidal  axis  from  the  center  of  curva- 
ture is  changed  from  R  to  p.  It  should  be  noted  that  y,  R,  and  p  at  any 
section  are  measured  from  the  centroidal  axis,  and  not  from  the  neutral 
axis. 

Let  ds  denote  the  length  of  the  centroidal  fiber  00'  and  e0  the 
unit  strain  of  this  fiber.    Then 


e0 


OiOx' 
00, 


or   0i0(  =  e0ds  =  eQRdO 


e  = 


PiPi' 
PPi 


PXH  +  HPX'  0i0(  +  HP( 


PPi 


PPi 
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But,  as  found  above,  OiO/  =  e0Rdd; 
further  HP/  =  Oi'H  Add  =  yAdd; 
also  PPi  =  (R  +  y)  dd 

Therefore  .  =    ***  +  »*» 


eQR  +  y 


(R  +  y)  dd 

Me 


de 


R  +  y 

For  convenience  let  the  angular  unit  strain  -=—  be  denoted  by  w. 

dd 

By  adding  and  subtracting  e0yi  the  expression  for  e  may  be  transformed 

to  the  following  equation : 

y 
e  =  e0  +  (w  -  e0)  ~]f~+~y  (2) 


From  equation  (1)  then 


[■ 


=  Ee  ■---  E  I  e0   h  («  -  c0) 


ft +  2/ J 


(3) 


Since  there  are  three  unknown  quantities  (<S,  w,  and  e0)  in  this  equation, 
two  additional  equations  must  be  found.  These  are  obtained  from  the 
conditions  of  equilibrium:  namely,  (1)  that  the  algebraic  sum  of  the 
external  forces  about  any  axis  in  the  cross-section  is  equal  to  the  mo- 
ment of  the  stresses  on  the  section  about  the  same  axis,  that  is,  the 

resisting  moment  is  equal  to  the  bending  moment,  or    /  syda  =  M; 

and  (2)  that  the  algebraic  sum  of  all  the  normal  stresses  on  the  cross-section 

is  equal  to  zero,  or  /  sda  —  0. 


Thus 

y 


M 


=   I  ysda  =    I  Ey   \  e0  +  (w  -  e0) 


da         (4) 


R+yJ 

I  sda  =   I  E\  e0  +  (<a  -  eo)  J^TT  \  da  =  0  (5) 

If  the  modulus  of  elasticity  E  has  the  same  value  for  all  the  fibers, 
it  is  a  constant  in  the  equations,  and  hence  equations  (4)  and  (5)  may  be 
written 


/ 
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M  =  E  T  e0   fyda  +  («  -  e0)  J  ~^  da  (6) 

and  e0  I  da=  —  (a-  e0)  I  —~  da  (7) 

But  /  da  —  a,  and  since  y  is  measured  from  the  centroidal  axis, 
yda  =  0.     Further,  for  convenience,  let 


/ 


- —  da  =  —Za 

R  +  y 


in  which  Z  is  a  property  of  the  area  somewhat  similar  to  the  moment  of 
inertia  in  the  straight-beam  formula.  The  integral  in  equation  (6)  may 
also  be  expressed  in  terms  of  Z  as  follows : 

f  v2  f  (        Ry  \  C   y 

I   77", —  da  =   I    [  y  —  77^ —  )  da  =  —  R  I   777 —  da  =  ZaR 

J  R  +  y        J    \       R  +  y  J  J  R  +  y 

Equations  (6)  and  (7)  now  become 

M  =  E  (w  -  e0)  ZaR  (9) 

e0  =  («  -  c0)  Z  (10) 

J_  /M       M  \ 

Substituting  these  values  in  equation  (3) 


M  M  1  1     (M  .    M\    (11) 

Hence  (w  -  c0)  =  "^^  •  ^o  =  ^  and  co  =  —  (  - 


S~    aR 


(12) 


in  which  S  is  the  unit  stress  at  a  point  at  the  distance  y  from  the  cen- 
troidal axis  of  the  section  (y  becomes  c  at  the  extreme  fibers) ;  M  is  the 
bending  moment  at  the  given  section ;  R  is  the  radius  of  curvature  of  the 
centroidal  axis  of  the  section  of  the  unstressed  beam;  a  is  the  area  of  the 
cross-section;  and  Z  is  a  property  of  the  cross-section  defined  above. 

The  bending  moment  M  is  positive  when  it  increases  curvature  and 
negative  when  it  decreases  curvature;  y  is  positive  when  measured  to- 
ward the  convex  side  of  the  beam  and  negative  when  measured  toward 
the  concave  side,  that  is,  toward  the  center  of  curvature. 
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Location  of  Neutral  Axis 


According  to  equation  (12)  the  unit  stress  on  any  fiber  at  the  dis- 
tance y  from  the  centroidal  axis,  due  to  the  bending  moment  M,  is 


S  = 


M. 

aR 


1  + 


Z\R  +  y) 


(13) 


The  stress  at  the  neutral  axis  is  zero,  and  if  y0  denotes  the  distance  of 
the  neutral  axis  from  the  centroidal  axis,  the  value  of  y0  may  be  found 
by  equating  S  in  the  above  equation  to  zero. 


Thus 


Z  \R  +  yQJ 


Vo  =   ~ 


ZR 
Z+l 


(14) 


The  negative  sign  means  that  y0  is  measured  from  the  centroidal 
axis  toward  the  center  of  curvature.  It  is  important  to  note  that  equa- 
tion (14)  applies  to  a  beam  that  is  subjected  to  bending  only. 

Axial  and  Bending  Loads  Combined 

If  a  curved  beam  is  subjected  to  a  normal  load  P  acting  through 
the  centroid  of  the  area  of  the  section  (axial  load)  in  addition  to  a  bend- 
ing moment,  the  resulting  unit  stress  is  the  algebraic  sum  of  the  unit 
stresses  caused  by  the  axial  and  the  bending  loads.  Thus  the  maximum 
unit  stress  is 


P        M 

s=  ~  ±  ~B 
a         aR 


(15) 


the  signs  of  M  and  y  are  determined  as  stated  above,  and  the  sign  of  P 
is  positive  when  it  produces  tensile  and  negative  when  it  produces  com- 
pressive stress. 
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